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Abstract. Caustics are high-density structures that form in collisionless media. Under self-gravity, cold dark 
matter flows focus onto caustics which are yet to be resolved in numerical simulations and observed in the real 
world. If detected, caustics would provide strong evidence for dark matter and would rule out alternative models 
such as those with modified dynamics. Here, we demonstrate how they might be observed in weak lensing data. We 
evaluate the shear distortion and show that its radial profile is marked by a characteristic sawtooth pattern due to 
the caustics in dark matter haloes that form by selfsimilar accretion. We discuss the observational complications, 
mainly due to the poor knowledge of the virial radii of the haloes and demonstrate that a superposition of about 
200 cluster-size haloes would give a signal-to-noise ratio which is sufficiently large for the detection of caustics 
with ground-based observations. This number is reduced to 60 for space-based observations. These bounds can 
be easily achieved by the ongoing wide field optical surveys such as CFHTLS and the future space-based projects 
SNAP and DUNE which have to be accompanied by an X-ray follow-up of the selected clusters for a precise 
determination of their virial radius. 
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, ^ 1. Introduction 



The nature of dark matter, which constitutes about 30% 
of the mass of the Universe, remains largely unknown. 
Results from cosmic microwave background explorers and 
large-scale galaxy surveys suggest that dark matter is 
cold with little velocity dispersion (e.g. Spergel et al. 
2003; Tegmark et al. 2004). If so, then its evolution 
is mainly governed by its self-gravity and expressed by 
Jeans- Vlassov-Poisson equations. The collisionless nature 
of dark matter predicts the formation of multi-stream re- 
gions bounded by very high density manifolds known as 
caustics. 

In a cold dark matter Universe, caustics can form on 
large scales of many megaparsecs, manifesting in the fil- 
amentary structure of galaxy distribution (Shandarin & 
Zeldovich 1989) and also at smaller scales of a few par- 
sees or kiloparsecs in dark matter haloes. This is well il- 
lustrated in some pionering numerical works (Melott & 
Shandarin 1989). Due to their abundance, the rich obser- 
vational and numerical data and their high density con- 
trasts, haloes are likely areas for the caustics. 
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Analytic models for the formation and evolution of 
dark matter haloes are still rare and most works are based 
on the selfsimilar accretion model (Gott 1975; Gunn 1977; 
Fillmore & Goldreich 1984; Bertschinger 1985). In this 
model, first proposed to explain the rotation curve of 
galaxies, haloes form by temporally self-similar collapse 
of dark matter shells onto an initially over-dense per- 
turbation. Dark matter shells initially expand until they 
reach their turnaround radii where they separate from the 
background expansion and collapse. After collapse they re- 
bounce and collapse again and the density profile settles 
asymptotically into a power-law which is convolved with 
singular spikes, namely with caustics. 

In the following we shall refer to such caustics as " outer 
caustics". They are suggestive of the sharp stellar shells 
observed around giant ellipticals which can arise in the 
merger of galaxies (Malin & Carter 1980; Quinn 1984; 
Fort et al. 1986). The main observational difference be- 
tween the merger and selfsimilar spherical infall is that 
the former predicts that the caustics are interleaved with 
the caustic radii alternating on opposite sides of the galaxy 
and the latter predicts concentric spherical shells. 

The spherically symmetric model has been extended to 
consider infalling matter with angular momentum and cal- 
culated the properties of an additional kind of "inner caus- 
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tics" with torus-like topology (Sikivie 1998, 1999). Such 
accretion with angular momentum is more relevant for 
galaxy-size haloes than for clusters of galaxies dominated 
by radial infall. 

Here, we focus on the outer caustics of cluster-size 
haloes and argue that they can be reasonably approxi- 
mated by selfsimilar infall models. In the original version 
of this model, dark matter is absolutely cold, i.e. with zero 
velocity dispersion, and caustics are infinitely thin concen- 
tric spherical shells with diverging densities. However, re- 
alistically, dark matter has a small velocity dispersion and 
these shells have finite thicknesses. The thicknesses of the 
caustics would however remain very small (due to the cold- 
ness of dark matter) and thus they would contain very lit- 
tle mass in spite of their significant density. Various char- 
acteristics of the caustics such as their density profile, their 
thickness and their approximate maximum density for a 
low velocity dispersion dark matter medium have been re- 
cently evaluated (Mohayaee & Shandarin 2006, hereafter 
MS05). 

Detection of dark matter caustics remains a challeng- 
ing problem both for 3-dimcnsional numerical simulations 
and for observations. Dark matter annihilation in the 
caustics has already been studied (See e.g. Sikivie ct al. 
1997; Hogan 2001; Mohayaee & Shandarin 2006; Pieri & 
Branchini 2005). Since the flux of the annihilation prod- 
ucts, e.g. 7-rays or antiproton flux, depends on the square 
of the local density 1 , caustics with their sharply-peaked 
densities would be the likely places for significant emissiv- 
ity Such kind of observations are promising for the even- 
tual detection of dark matter in the caustics. However, 
the major problem with dark matter detection through 
gamma-ray emission is the severe background contami- 
nation of the signal and the low signal-to-noise ratio of 
present day observations. 

It has been shown that rotation curves of galaxies 
might be sensitive to the presence of inner caustics and 
claimed a marginal detection based on an ensemble av- 
erage over 32 rotation curves (Kinney & Sikivie 2000). 
Inner caustic rings are located in the plane of the disk 
and are likely to modify rotation curves more efficiently 
than outer caustics with spherical symmetry. Moreover, 
rotation curves only probe the dark matter potential of 
galactic haloes where a large amount of tracers (gas or 
stars) is available. Therefore, outer caustics which are lo- 
cated well beyond the observable tracers cannot leave a 
detectable imprint on rotation curves. 

Gravitational lensing provides a promising alternative 
tool. Since lensing probes the projected density profile 
with no regards to the nature or dynamical state of the 
deflecting mass, it should be sensitive to the caustics. The 
lensing properties of dark matter caustics and particularly 
their efficiency in magnifying and/or producing multiple 
images of background sources has already been investi- 
gated (Charmousis et al. 2003). 

1 Sec e.g. Donato et al. (2004) for the antiproton flux and 
Bertone & Merritt (2005) for a general recent review. 



Inner caustics may be dense enough to produce substan- 
tial magnification and small separation multiple images as 
in micro-lensing by compact objects. High magnification 
events due to such caustics may explain the anomaly in 
flux ratios observed in multiply-imaged quasars which are 
hardly reproduced by a smooth halo or even subhalos (see 
e.g Dalai & Kochanck 2002; Kochanek & Dalai 2004). 
Conversely, it has been found that outer caustics are in- 
efficient in magnifying distant sources and would yield 
at most a few percent net magnification or shear 
(Charmousis et al. 2003). As we shall detail below, the 
selfsimilarity of dark matter accretion implies that outer 
caustics occur at the same radius provided physical radii 
are properly rescaled by the halo mass. Hence the tiny 
lensing signal of caustics could show up statistically by 
averaging over many rescaled haloes. Such a statistical 
approch is much more complicate for the lensing-based 
detection of inner caustic rings because of the random 
orientation of rings (or angular momenta) . 

In this work, we consider the lensing properties of the 
outer caustics only and we propose the weak-lensing ef- 
fect as a potential way to detect caustics. We demonstrate 
that the caustics will produce sharp variations in the pro- 
jected surface mass density around haloes. Depending on 
the height and width of caustics, gravitationally-distorted 
background galaxies will experience local variations of 
shear. 

If the aforementioned universal property of haloes is 
fullfilled, then the rapid progress in X-ray and lensing ob- 
servations of cluster of galaxies may offer one possibility 
to observe dark matter caustics. Deep and wide surveys 
such as CFHTLS will cover fields of view of a few hun- 
dreed square degrees and will provide us with a useful 
material for the detection of dark matter caustics through 
the capability of mass to coherently stretch the image of 
the background galaxies. They will provide us with a large 
enough number of galaxy clusters 2 to achieve the required 
level of signal-to-noise ratio. Wide field spatial surveys 
will be even more powerful for the investigation of the 
lensing properties of dark matter haloes and their associ- 
ated caustics. An important requirement for the detection 
of caustics is the measurement of virial radii of the clus- 
ters which can be provided by X-rays (e.g. Arnaud et al. 
2005). 

Throughout this paper, we assume an Einstein-de 
Sitter Universe but our results should be qualitatively sim- 
ilar in a concordance ACDM model. The role of dark en- 
ergy becomes important at low redshifts (<~ 0.2) which 
we expect to occur well after the formation of the typical 
dark matter haloes we consider here. Furthermore, once a 
particle turns around and collapses, it separates from the 
background expansion and its subsequent motion should 
not be affected by the A term. 

The paper is organized as follows. We review the three- 
dimensional and projected properties of self-similar haloes 
and caustics in Sect. 2. We derive the lensing signal for a 

2 ~ 5 per square degree, (Hennawi & Spergel 2005) 



Raphael Gavazzi et al.: Probing dark matter caustics with weak lensing 



3 



single halo, compare it to the noise level of fiducial ob- 
servations and estimate the number of haloes required to 
achieve a significant signal-to- noise ratio in Sect. 3, where 
we also examine the ability of weak lensing to constrain 
the velocity dispersion of dark matter particles. We sum- 
marize, discuss the prospectives for future works and con- 
clude in Sect. 4. 

2. Dark matter caustics 

2.1. Tri-dimensional key equations 

In an Einstein-de Sitter Universe a spherical overdensity 
expands and then turns around to collapse. After collapse 
and at late times, the fluid motion becomes selfsimilar: 
its form remains unchanged when its length is re-scaled 
in terms of the radius of the shell that is currently at 
the "turn around" and is falling onto the galaxy for the 
first time. Physically, selfsimilarity arises because grav- 
ity is scale-free and because mass shells outside the initial 
overdensity are also bound and turn around at successively 
later times. Self-similar solutions give power-law density 
profiles whose exact scaling properties depend on the cen- 
tral boundary conditions and on whether the fluid is col- 
lisionless or collisional (Gott 1975; Gunn 1977; Fillmore 
& Goldreich 1984; Bertschinger 1985). The density profile 
obeys a power-law on the scale of the halo which provides 
an explanation of the flattening of the rotation curves of 
the galaxies. However, on smaller scales the density pro- 
file contains many spikes (i.e. caustics) of infinite density. 
The position and the time of formation of these caustics 
are among the many properties that have been studied in 
the selfsimilar infall model (Bertschinger 1985). 

In the presence of a small velocity dispersion the max- 
imum density and thickness of the caustic shells and their 
density profiles have been evaluated in the framework of 
a selfsimilar collapse model (MS05). 

The global halo density profile, asymptotically reached 
in this process, is well-approximated by 



Phalo(A) 2.8 A- 9 / 4 



p H (1 + A3/4)2 ' 



(1) 



where A = r/r ta with r ta the present turnaround radius 
of the halo and pn = 3H 2 /8ttG the background density. 
The turnaround radius can be easily computed using the 
virial radius. Within the virial radius, r v ; r = ^200, the 
mean density, p v - 1T , is, by definition, 200 times the back- 
ground density. Thus, using the density profile (1), we 
obtain the following relationship between the turnaround 
and the virial radii: r ta ~ 4r v i r . 

For a perfectly cold dark matter medium, the density 
profile close to a caustic at Afe is (Bertschinger 1985) 



Po (A) 



Gk 



PH 



VAfc - A ' 



a = 0. 



with 



Gk = 



4 V /T 2AI 



-2? fc /3 

"at 



(2) 
(3) 



where the values of the various quantities Afe, [and 
Afc which will appear in the coming expression (6)] are 
given in Table 2.1 [see MS05 for a detailed description of 
these parameters]. 

When the temperature of particles is not strictly zero 
and the velocity of particles is distributed according to the 
distribution function f(v), caustic positions are shifted by 
a small value SX and the caustic density is modified as 
(MS05): 



Pa {X) = Jdvp [X- SX(v)} f(v) ; a ± 



(4) 



In this work, we choose a top-hat velocity distribution 
function (MS05) . Then, the density close to the fc-th caus- 
tic in the halo is given by 



MA) 
Ph 



Gk 
Ak~ 



X for A < A fe , 

for X- < X < A+ , 
for A > A£ , 



(5) 

where A fe = Afe — Afe and A^ = Afe + Afe and the thickness 
of the fc-th caustic, Afe, in nondimensional coordinate is 
given by 



nta 



Afe = Arfe = i^j_ e 5W9 | A, | ta(t) (6) 



and nta is the initial turnaround radius and a(t) is the 
value of the velocity dispersion of dark matter particles at 
time t which is that at decoupling re-scaled by the ex- 
pansion factor 3 . The standard spherical collapse model 
yields a relation r;t a = 2r v i r for a constant overdensity 
[see e.g. chapter 5.10 of (Padmanabhan 2002)], However, 
the real value of the initial turnaround radius would be 
lower due to the continuous accretion by the halo. Thus 
the physical thickness of the caustic, Ar fe , depends only 
on its position in the halo and the nature of dark matter. 

The maximum density at each caustic position can also 
be approximated by the expression (see Table 2.1 and also 
MS05) 



pmax — 



(37r)i/3^AT V 1 ^) 



2Gfce- 5 «*/ 18 njl 



PH 



(7) 



Expressions (5), (6) and (7) provide us with a sufficient 
basis for the evaluation of magnification and shear due to 
dark matter caustics. Although in the rest of this work we 
directly integrate expression (5) and never use (7), values 
for the latter are given in Table 2.1 to demonstrate the 
density contrast of each caustic with respect to its host 
halo. 

Although the selfsimilar model might seem naive, 
it provides a good approximation for outer caustics in 
galaxy-cluster haloes which are not significantly disrupted 
by merger and substructures and to a good approximation 

3 Hereafter, velocity dispersion is given at the present time, 
z = 0. For instance, neutralinos have a ~ 0.03 cm/s and axions 
have a ~ 10~ 7 cm/s. 
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(d 2 \/df) k 


Afe 


Ar k (pc) 


/ / / i \ — 1/2 

Pmax/pJ//(V r vir) pC ' 


phalo/p// 


1 


0.985 


0.368 


-5.86 


-0.0704 


5.6 10" 4 


233 


12 


2 


1.46 


0.237 


-11.2 


-0.0254 


2.6 10~ 4 


431 


39 


3 


1.76 


0.179 


-16.7 


-0.0135 


1.7 10~ 4 


640 


82 


4 


1.98 


0.146 


-22.3 


-0.00854 


1.2 10~ 4 


850 


139 


5 


2.16 


0.124 


-28.0 


-0.00591 


9.1 10~ 5 


1067 


210 


6 


2.31 


0.108 


-33.9 


-0.00437 


7.3 10~ 5 


1290 


297 


7 


2.43 


0.096 


-39.8 


-0.00337 


6.0 10~ 5 


1532 


398 


8 


2.55 


0.087 


-45.7 


-0.00266 


5.1 10~ 5 


1749 


513 


9 


2.64 


0.079 


-51.7 


-0.00221 


4.4 10~ 5 


2010 


640 


10 


2.73 


0.073 


-57.8 


-0.00182 


3.8 10~ 5 


2253 


785 



Table 1. The non-dimensional times (£fc), positions (Afc), thicknesses (Ark), maximum densities (pmax) of the first ten caustics, 
the halo density itself (phalo) at the positions of these caustics and other parameters used in expressions (3), (5), (6) and (7). 
The present velocity dispersion is that for neutralinos which is about 0.03 cm/s. It is worth emphasizing that the thicknesses of 
the caustics do not depend on the halo parameters but only on the velocity dispersion of dark matter [see expression (6)]. It is 
instructive to compare the last two columns, which show that for a cluster-size halo with r v i r ~ 1 Mpc, the caustic density can 
be many orders of magnitude higher than the halo density itself. 



are spherical. Furthermore, fluctuations caused by large 
scale structure would already be averaged out in the sta- 
tistical evaluation of the shear. 



2.2. Projected densities 

Since we are concerned with the lensing properties of caus- 
tics, we have to calculate the projected density profile for 
the caustics and for the halo. The Abel integral relates the 
3-dimensional density (p) and the 2-dimensional density 
profiles (£) by 



E(A) = 2r taj 



p(X')X'd\' 
VA' 2 - A 2 



(8) 



We numerically integrate the above expression, for the 
density profiles (2) and (5). The halo surface mass density 
can be evaluated analytically if we neglect the (1 + A 3 / 4 ) 2 
term in equation (1), i.e. at small scales A <C 1, yielding 
the approximate projected halo profile 



£ halo (A) ~ r ta 7.56\-^ 4 p H ■ 



(9) 



Let us define the mean projected density enclosed by 
the cylinder of radius A: 



v(A) = A / £(A')A'dA'. 



(10) 



Thus, under the hypothesis of small A, the mean projected 
density for the halo component is £haio(A) = §S(A). 

Figure 1 shows the 3D (left) and 2D (right) density 
profiles for the halo and the caustics. We consider the case 
of a perfectly cold dark matter with very peaked caustics 
and the case in which they are smoothed out by a finite 
velocity dispersion a <~ 60kms _1 . This value is very high 
for most cold dark matter models but its extremity illus- 
trates well the peak dilution and shows that caustics sur- 
vive a large thermal softening. When considering the pro- 
jected density, instead of spikes we rather see stairs that 



come from the analytic 1/V r fc — r singularity of caustics. 
Consequently even with singular caustics, the projected 
density profile is not peaked. The implications for lensing 
are discussed in the next section. 

3. Weak lensing 

3.1. General equations 

The fundamental quantity for gravitational lensing is the 
lens potential tj){6) at the angular position 6 which is re- 
lated to the surface mass density £(0) projected into the 
lens plane through : 



4G A A 



J & 2 9'Y>{0')\n\0-0'\ , (11 



where A, A and As are angular distances to the lens, 
to the source and between the lens and the source respec- 
tively. The deflection angle a = V-0 relates a point in 
the source plane f3 to its imagc(s) in the image plane 9 
through the lens equation /3 = 6—a.(6). The local relation 
between /3 and is the Jacobian matrix Aij = dfii/dOj : 



1 - K - 71 

-72 



-72 

1 - K + 71 



(12) 



with the convergence k(6) = £(0)/£ cr it directly related 
to the surface mass density via the critical density 



Ac 



4ttG A A 



(13) 



and the 2-component shear 7 = 71 + 272 in complex 
notation. The convergence satisfies the Poisson equation 
Aip = 2k. In the weak lensing regime (7 <C 1), an ellipti- 
cal object in the source plane with complex ellipticity e s is 
mapped into an elliptical image with a different ellipticity 
e = e s + 7. We refer the reader to the reviews of Mellier 
(1999) and Bartelmann & Schneider (2001) for detailed 
accounts of weak lensing. 
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Fig. 1. Left panel: In units of the background density, pH, we show the 3D density profile for the halo (black), for the caustics 
in a perfectly cold medium (red) and for the caustics smoothed out by a warm/hot dark matter with velocity dispersion 
a = 60kms _1 ( green). Red spikes are all singular with infinite density and limited here to finite values because of finite 
resolution. Right panel: The same colour-coding for the projected 2D density profiles in units of puTta- When projected, 
caustic peaks are smoother and look like a flight of stairs. Therefore, we expect low lensing magnifications close to the caustics 
(see Section 3). 



For a circularly symmetric lens, 7 is oriented tangen- 
tially to the lens center and its amplitude at radius r is 
7(r) = (S(r)— X(r))/£ cr i t . Since sources are randomly ori- 
ented, the tangential component of the observed galaxies 
is an unbiased estimator of 7. When averaging the esti- 
mate of 7 within an aperture of solid angle S7, containing 
N — fin galaxies (n is the number density of sources), 
the noise dispersion of 7 is 07 = -j= where a e ~ 0.3 is 
the intrinsic dispersion of source ellipticitics (along one 
component). 

3.2. Shear measurement 

In order to be consistent with our calculations of Sec. 2.2, 
we define a pseudo-shear: T(A) = (S — £)/ ' pu^ta and the 
corresponding noise level rjv = S C ritC 7 //0/frt a . For an EdS 
cosmology, and considering an annulus of inner and outer 
radii Ai and A2 respectively, it is straightforward to write 
Tat in this useful form: 

r,(A 1 ,A 2 )=2.16^f^) 2 (l + , i) 3 x 

L>U \ T ta J 

/30arcmin -2 / a e \ 1 ,„ 

V— 5— (os)^^- (14) 

This expression seems to indicate that the noise level will 
be lower for nearby haloes. However, it hides the fact 
that low redshift haloes require a very wide sky cover- 
age for the outermost caustics (~ r ta /3) to fit the field of 
view of the observation. Consequently, intermediate red- 
shift haloes (z ~ 0.2 — 0.5) are the most interesting targets. 
In addition, nearby (and thus large angular scale) clusters 
suffer from noise due to large-scale structure (LSS) fluc- 
tuation integrated along the line-of-sight and unrelated to 
the halo we are considering (Hoekstra 2003). For scales 



< 15 arcmin, the smearing of the shear profile by LSS is 
minimised. 

Figure 2 shows T as a function of A for the same val- 
ues of the thermal velocity dispersions a = Okms -1 and 
a = 60kms~ . Comparing these curves, one can see that 
the sawtooth patterns due to caustics survive significantly 
high temperatures. Next, we consider the noise level for a 
fiducial halo at redshift z\ — 0.3 and a turnaround radius 
r ta = 5 Mpc which is a typical value for clusters (upper 
green "binned" curve). With a single halo the detection 
of caustics is impossible. If we are able to stack the signal 
from a few tens of clusters (upper blue-binned curve) , the 
noise level will be low enough to be sensitive to caustics as 
a whole. However, we would detect a smooth contribution 
which is indistinguishable from the halo itself. 

Sawtooth patterns cannot be confused with the effect 
of substructures since the contribution of the latter would 
be averaged out over the azimuthal angle (most outer 
caustics have spherical symmetry) and once rescaled, caus- 
tics always appear at the same place within haloes. This 
is not the case for substructures which can appear at any 
radius inside the host halo. 

Consequently, the right way to probe the existence of 
caustics is to measure the T signal in excess/default rel- 
ative to the extreme value of a, e.g. a > 300 km s 1 as 
taken here. For this purpose, one needs > 100 stacked 
clusters. With the corresponding noise level, it would be 
possible to test the thermal smoothing of caustics and put 
constraints on a. However, the sensitivity is poor and only 
upper limits can be put on a with a realistic number of 
haloes. For instance with TV = 100 (resp. 250) clusters, we 
could achieve a limit a < 170 km s -1 (resp. 40kms~ 1 ) at 
a 95.4% confidence level. 

When considering galaxies instead of clusters with 
turnaround radii about 10 times smaller, the number of 
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Fig. 2. r(A) contribution of caustics for two values of a: cold 
medium a — (black curve) and warm medium a ~ 60km/s 
(red curve). The green (resp. blue) "binned" curve is the noise 
level for one (resp. 100 stacked) halo(es). 



haloes required to achieve the same detection level is 
10 4 times higher. So using a few 10 5 galaxies between 
z <~ 0.1 — 0.5 would yield the same results. 

This required level of signal can easily be achieved 
with a wide and moderately deep survey like the ongo- 
ing CFHTLS. Typically, a square-degree field of view will 
contain a few such cluster-size haloes with r ta ~ 5 Mpc, 
a few thousand elliptical galaxies with r ta ~ 1 Mpc and 
tens of thousands of spiral galaxies with r ta ~ 500 kpc. 
The total coverage of the CFHTLS wide survey is 170 
square degree and will contain a large enough number of 
clusters/galaxies. Furthermore the wide fields of view are 
well-suited to measure shear up to the outermost caustic 
of clusters (1 Mpc = 5' at z = 0.3 and 4' at z = 0.5). 

Spatial observations provide an improvement on shear 
measurements: (i) the intrinsic dispersion in source ellip- 
ticities is lowered a e ~ 0.2 and (ii) the density of sources 
is increased n ~ 50 axcmin . Hence the total number 
of haloes required to achieve the same detection level can 
be lowered by a factor of three. Future wide spatial sur- 
veys like SNAP or DUNE will provide the required sky 
coverage. 

3.3. Halo-stacking issues 

So far, we have considered the observational difficulties 
encountered by the dispersion in intrinsinc source elliptic- 
ities. However, practical complications such as the diffi- 
cult signal stacking process should also be overcome be- 
fore a detection of the caustics can be achieved. To do so, 
high-precision measurements of the location of the cen- 
ter and the turnaround radius of each halo are required. 
Otherwise, imperfect alignment/rescaling would tend to 
blur the caustics spikes and reduce the sensitivity. The 
center of a well-relaxed cluster coincides within one arc- 
second with the center of the halo. Hence the only scaling 



factor is the turnaround radius that can be related to the 
virial radius [see the sentence following expression (1)]. 




Fig. 3. Difference between the P(A) contribution of cold caus- 
tics and the smooth component. We illustrate the blurring ef- 
fect of an imperfect knowledge of the turnaround radius of 
stacked haloes. The solid black line, dashed green line and dot- 
dashed blue line correspond to zero, 3% and 10% uncertainties, 
respectively. As in Fig. 2, the binned curve represents the noise 
level achieved with 200 stacked haloes. The convolution effect 
of error in the assumed/measured value of r ta is important 
and a significant detection of caustics requires well-determined 
turnaround radii (< a few percent relative accuracy). The very 
small scale oscillations in the plot are numerical artifacts and 
ideally the only sawtooth patterns are those due to caustics. 

Figure 3 shows the blurring of caustics due to an im- 
perfect knowledge of the turnaround radius of each stacked 
halo. We consider a dispersion around the true value of 3% 
and 10%. To properly detect caustics, one needs a precise 
estimate of r ta . Provided dark matter is cold enough, a < 
a few kms -1 , which is a realistic prescription, a reason- 
able number of stacked haloes, ~ 200, can overcome the 
loss of a few percents of relative precision caused by the 
errors in the determination of the turnaround radius. 

Consequently the number of haloes necessary for de- 
tection needs to be increased in order to achieve the 
necessary level of signal-to- noise ratio. Moreover, addi- 
tional external information for the determination of r ta 
such as X-rays or dynamical observations in addition to 
weak lensing data are indispensable for a convincing de- 
tection of caustics. The present-day state-of-the-art X-ray 
estimates for the scale radius is < 10% for nearby clus- 
ters (z < 0.2) (Arnaud et al. 2005; Pointecouteau et al. 
2005) and it seems that a similar precision can be reached 
at higher redshifts (0.4 < z < 0.7) (Kotov & Vikhlinin 
2005). In addition shear, which is used for caustic detec- 
tion, also provides constraints on the halo density profile. 
The virial radius of some clusters presenting strong and 
weak lensing features can be measured with good accu- 
racy (Ar 20 o/r2oo ^ 3%) (Broadhurst et al. 2005; Gavazzi 
2005). Consequently, future large cluster samples with X- 
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rays and lensing data of present-day quality will provide 
us with the necessary precision to probe dark matter caus- 
tics. 

4. Discussion & conclusion 

Although the Liouville theorem claims that singularities 
will survive, they are likely to develop a complex topology 
in the course of evolution of structures under self-gravity. 
It is not clear to what degree the merging processes and 
substructures will smear out the caustics or complicate 
their geometries. Here, we have considered the outer caus- 
tics of cluster-size haloes which are expected to have suf- 
fered far less from mergers and due to their relatively large 
separations are unlikely to have been washed away by the 
dispersive-like effect of the substructures. These caustics 
are the ones that contribute most to the lensing signal 
for the following reasons. Their amplitude relative to the 
background smooth halo component is more important as 
compared to the inner caustics. The inner caustics smear 
most from thermalisation and imperfect stacking and fi- 
nally, the noise level increases toward the center of haloes 
thus giving more weight to the outer caustics. Hence, we 
have focused on caustics of cluster-size haloes and have 
argued that they can be reasonably approximated by a 
sclfsimilar spherical accretion model, though the triaxial- 
ity of haloes is well established in numerical simulations 
(Jing & Suto 2002). However, caustic patterns should ex- 
ist in triaxial matter distributions. The singularities will 
have the same elliptical symmetry and may be properly 
stacked from one halo to another by choosing a subsam- 
ple of apparently circular projected haloes or by using the 
shear azimuthal variations to constrain the halo ellipticity. 
We expect this effect to be comparable to the uncertainty 
in the halo scaling radius (cither turnaround or virial). 

We have shown that the existence of dark matter caus- 
tics could be probed by properly stacking the weak lens- 
ing signal of a reasonable number of haloes. The main 
observational limitation is perhaps the precise estimation 
of the turnaround radius, r ta , of superimposed haloes but 
we have shown that the loss of a few percent relative ac- 
curacy in the determination of r ta (or asphcricity) can be 
compensated for by stacking about 200 haloes. 

Although the sensitivity is low, a detection of caus- 
tics provides an upper bound for the temperature of dark 
matter, thus excluding hot dark matter models. The sen- 
sitivity is not sufficient to distinguish between various cold 
dark matter candidates (like axions or ncutralinos) since 
for most of the corresponding velocity dispersions, the 
shear signal would be similar. 

However, a detection of caustics would be a strong argu- 
ment for the existence of cold dark matter since alternative 
models like MOND could not explain such density singu- 
larities and at most could serve in place of a smooth halo 
(namely provide an equivalent effective gravitational po- 
tential). 

Putting constraints on the velocity dispersion of dark mat- 
ter is a challenging topic in modern cosmology since it 
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offers the possibility to pin-down an actual physical pa- 
rameter of dark matter. In this paper, we investigated the 
possibility of such a measurement with the weak lensing 
effect of dark matter caustics. The implicit observational 
hypothesis is that we can have a sclfsimilar geometrical 
description of the caustic shell distribution which depends 
on a single characteristic scaling parameter: i.e. the virial 
radius. 

Wide field surveys such as the ongoing CFHTLS ac- 
companied by X-ray observations can provide the required 
statistics for a successful detection of caustics. The num- 
ber of haloes required to be superimposed will be lowered 
by a further factor of 3 for future space-based experiments 
like SNAP or DUNE. 

Here, we have used the first and most common caus- 
tic singularity, that for which the density profile falls 
with inverse square-root of the distance from the caus- 
tic. Caustics of higher-order singularities can also appear 
in collisionlcss media and have already been classified 
(Arnold 1986). It remains a challenging task to generalize 
our work to haloes with non-vanishing eccentricity where 
higher-rank caustics would occur and to examine if they 
can modify the magnification properties of the lensed im- 
ages and account for anomalous image flux ratios. 
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